BLOCK-DIAGONAL SEMIDEFINITE PROGRAMMING 
HIERARCHIES FOR 0/1 PROGRAMMING 
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Abstract. Lovasz and Schrijver, and later Lasserre, proposed hierarchies 
of semidefinite programming relaxations for general 0/1 linear programming 
problems. In this paper these two constructions are revisited and two new, 
block-diagonal hierarchies are proposed. They have the advantage of be- 
ing computationally less costly while being at least as strong as the Lovasz- 
Schrijver hierarchy. Our construction is applied to the stable set problem and 
experimental results for Paley graphs are reported. 

1. Introduction 

A basic approach in combinatorial optimization consists of formulating the prob- 
lem at hand as a 0/1 linear programming problem, typically of the form 

max {c T x \Ax<b, x£ {0,1}™}, 

where c <E R", b € M' m and A € M mx ". Then the task is to find an efficiently 
computable outer approximation of the polytope P, defined as the convex hull of 
the 0/1 solutions to Ax < b. 

On the one hand, extensive research has been done for finding (partial) lin- 
ear inequality descriptions for many polyhedra arising from specific combinatorial 
optimization problems. On the other hand, researchers focused on developing gen- 
eral purpose methods for arbitrary 0/1 linear programming problems. Here let 
us mention the method of Gomory for generating cuts strengthening the initial 
linear relaxation {x £ R™ | Ax < b} of P and its various extensions for gener- 
ating strong cutting planes (see e.g. [H] HE]), the lift-and-project method [I], the 
reformulation- linearization technique [17j . the matrix-cut method of Lovasz and 
Schrijver |13) . and the sums of squares and moment method of Lasserre [8]. Some 
of these methods are compared in [9] ; see also [11] . A common feature of the meth- 
ods of Lovasz-Schrijver and of Lasserre is that they consider hierarchies involving 
semidefinite relaxations of P: Convex sets Qt (t — I, . . . ,n + 1) are constructed 
which can be described by semidefinite conditions and which form a hierarchy of 
increasingly stronger relaxations: 

{x e R n | Ax < b} 2 Qi 2 Qi 2 • . . 2 Qn+i = P. 

The two hierarchies are related; it is shown in [9] that the hierarchy of Lasserre 
refines the hierarchy of Lovasz-Schrijver. 



Date: August 19, 2008. 

1991 Mathematics Subject Classification. 90C22, 90C27. 

Key words and phrases. 0/1 linear programming, semidefinite programming, stable sets, Pay- 
ley graphs. 

This work was supported by the Netherlands Organization for Scientific Research under grant 
NWO 639.032.203. The third author was partially supported by the Deutsche Forschungsgemein- 
schaft (DFG) under grant SCHU 1503/4. 



2 



N. Gvozdenovic, M. Laurent, F. Vallcntin 



In this note we revisit these hierarchies and propose two new ones, which differ 
in the way of encoding the linear constraints defining the starting linear relaxation 
of P. Moreover one of them (introduced in Section 12. 4p can also be defined when 
the starting relaxation of P is an arbitrary convex body, as is the case for the 
Lovasz-Schrijver construction. The new hierarchies are nested between the Lasserre 
and Lovasz-Schrijver hierarchies, but they are less costly to compute. So they are 
especially well-suited for implementations. E.g., at given order t, the new hierarchy 
from Section HOI involves l/(t + l)!n t+1 + 0(n r ) variables compared to 2 t ~ 2 n t+1 + 
O(n') variables for the Lovasz-Schrijver hierarchy and to 0(n 2t ) variables for the 
Lasserre hierarchy. The new hierarchies can be seen as a variation of the Lasserre 
hierarchy, where one replaces a large matrix of order O(n') by smaller blocks of 
order n + 1 arising by block-diagonalizing suitably defined principal submatrices of 
the original large matrix. The motivation for considering block matrices is that it is 
computationally easier to solve a semidefinite program involving many small blocks 
rather than one large matrix. Most currently available interior-point algorithms for 
semidefinite programming are indeed designed to exploit block-diagonal matrices. 
While the hierarchy of Lovasz and Schrijver is originally defined recursively, we 
give an explicit description obtained by "unfolding" the recursion. In this way, the 
connection to the new hierarchies becomes transparent (see Section [5] for details). 

When applied to the stable set problem, our new construction gives a block- 
diagonal hierarchy whose first two steps were already used in the literature. The 
first order relaxation gives the Lovasz theta number and the second order one gives 
parameters considered in [U [10] for the stable set problem and in [6l [7] for the 
coloring problem. In these applications the computational advantage of the new 
hierarchy was of crucial importance. 

Contents of the paper. In Section [2] we first briefly introduce the constructions 
of Lovasz-Schrijver and of Lasserre. Then we give the new construction and show 
how to derive more compact formulations by block-diagonalization. In Section[3]we 
apply it to the stable set problem and in Section U we present some computational 
results illustrating the behavior of the new hierarchy for approximating the stability 
number of Paley graphs. 

Notation. Given a finite set V, we denote the collection of all subsets of V by 
V(V). Given a non-negative integer r, set P r (V) := {I € V(V) \ \I\ < r} and 
V =r (V) :— {I £ V(V) I |/| = r}. By we denote the empty set. Sometimes we 
identify V = i(V) with V, i.e., we write i instead of {£}. Furthermore, we sometimes 
write ij instead of {i,j} and ijk instead of {i,j, k}, etc. The standard unit vectors 
in R' Pl ( 1/ ) are denoted by eo, &% for i 6 V. 



2. Semidefinite programming hierarchies 

Suppose we are given a convex cone K contained in the homogenized unit cube 
{x e | < asj < x (i € V)}. Set 

P K := conv jx e {0, l} y | Q e K 
C K ■= K+ | ( M £ K with x e {0, l} 1 
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The general objective is to find the linear inequality description of the polytope Pk 
or, equivalently, of the cone Ck- In Section [2TT] we recall the construction of Lovasz- 
Schrijver which applies to any convex cone K. While the original construction 
is recursive we propose an explicit semidefinite programming reformulation. In 
Section 12.21 we recall the construction of Lasserre which applies to the case when 
K is represented by polynomial inequalities. Here we focus on polyhedral cones K 
of the form 

K = {x G R Pl(v) | af x > (I = 1, . . . , m)}, (1) 

where a\, . . . , a m G R^ 1 ^). In Section ^. 31 we introduce our new construction, which 
can be seen as a variation of the previous methods. We discuss two new hierarchies. 
The first one applies to polyhedral cones K as in {T]) and is more economical than the 
Lasserre hierarchy while still refining the Lovasz-Schrijver hierarchy. The second 
one applies to any convex cone K and can be seen as a non-recursive analogue of 
the Lovasz-Schrijver hierarchy having a more compact and explicit formulation. 

2.1. The Lovasz-Schrijver hierarchy. In this section we recall basic facts about 
the Lovasz-Schrijver hierarchy. For proofs and more details we refer to [13] . Set 

M+y ■= {Y G R^OTi^v) \ YhO,Y ii = Y oi (i G V)}, 

where "y 0" stands for "is positive semidefinite" . For a convex cone K C R^ 1 ( y ) 
define 

M+{K) := {Y G M+y | Ye, G K, Y(e - e<) G K (i G V)}, 
N + {K) := {Ye \ Y G M+(K)}. 

The t-th iterate of the Lovasz-Schrijver hierarchy is Nl_(K) := N + (N^ 1 (K)) for 
t > 1, where N\{K) := N+(K) and Nl(K) := K. It lies between K and C K and 
iV| +1 (A') C N\{K). We have N\(K) = C K for t = \V\. Moreover, for any fixed t, 
if one can optimize over K in polynomial time then the same holds for N+(K). 

In the following proposition we "unfold" the recursive definition of Ni_(K) and 
give an explicit semidefinite programming formulation. Its proof is straightforward 
and thus omitted. 

Proposition 2.1. A vector x G R^ 1 ^ lies in N^_(K) if and only if there exist a 
matrix Y G -M+y and matrices y? 1 '"'!"' G M+y , with s = 1, . . . , t— 1, i\, . . . , i s G 
V and <7i, . . . , <T a G {±1} satisfying the following conditions: 

(a) x = Ycq. 

(b) For all s = 1, . . . ,t — 1, ii, . . . , i s G V, and a%, . . . , o~ s -\ G {±1}-' 

v CTl,...,a-„_l _ - v rO-l,...,<7 8 _l, + l 

J ii,...,i s _i ~~ 1 i 1 ,,,,,i s L 0, 

where X-' 1 '"; 1 '" 1 = V for s = 1. 

(c) for aZZ ii, . . . ,ij G V and <ti, . . . , o t -\ G {±1}' 

cxr^e*, 
cxr^-ejG^, 

w/iere y CTl --- '*- 1 = y for t = 1. 
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The above formulation allows to estimate the cost of optimizing over N!i(K) 
in terms of n = \V\ and t. Set h(n,t) : £*~o(2n) s = 2 t ~V~ 1 + 0(n'- 2 ). The 
formulation involves (^)h(n,t) variables, h(n,t) matrices of order n + 1, and (2n)* 
conditions of type "a; G K" . Furthermore, it turns out that for 1 < r < s, the z r -th 
column of the matrix Y? 1 ' ' ' '! <Ts is identically zero whenever ay = — 1 and it is equal 
to the 0-th column whenever ay = 1. Thus, Y^ 1 '"'^ 3 can be assumed to have order 
n — s + 1 . 

2.2. The Lasserre hierarchy. In this section we recall some basic facts about 
the Lasserre construction, applied to the case when K is of the form {T]); for more 
information we refer to [8] and [9]. The Lasserre hierarchy involves moment ma- 
trices: A matrix whose rows and columns are indexed by a subset A of V{V) 
is said to be a moment matrix if the (/, J)-th entry depends only on the union 
I U J (for all I, J £ A). In this definition one may allow A to be a multiset, 
which corresponds to repeated rows and columns in the moment matrix. For a 
non-negative integer t and a vector y G R^ 2 *^, define the moment matrix of y of 
order t by M t (y) :— (yiuj)i.jev t (v)- F° r a se t T and a vector y G Bf p ' T * ) , we write 
M T {y) := (yiuj)i,Jev(T)l thus M n (y) = M v {y) if |V| = n. 

The following fact, observed in [131 [9], explains the relevance of moment matrices 
to 0/1 polyhedra: For x G R y define Q x :— (Ilie/ x i)iev(v)- Then we have for 

M v (y) t y e K+ {Cx | x € {0, 1} V } 

^vscv. Yl (- 1 ) |s ' vs| ys' > o. ( 2 ) 

In Lemma 12.21 we give an extension of this result. 

Next we explain how to encode the linear constraints afx > describing K . 
Given y e R p 2 t (V) an d a E « Vl( ' V \ define the vector ay G R^-iiv) by ( ay ^ ._ 
ooy/ + X^ev a i2//u{i} f° r ^ € V2t-i(V). For t > 1 we define the iterate of the 
Lasserre hierarchy by 

Q<(#) := {x G | 3y G R p2t(y) : y = x , Vi = Xi (i G V), 

Mt(j/)hO, M t _i(aij/) b (/ = l,...,m)}. '' >l 

It lies between K and Ck and Q* +1 (i4r) C Q f (K). The Lasserre hierarchy refines 
the Lovasz-Schrijver hierarchy, since we have Q t+l {K) C N + (Q t (K)) which implies 
Q t+1 (K) C iV|(iT) and Q n+1 (iT) = C K . The formulation © involves J2?=o (?) = 
0(n 2t ) variables, one matrix of order J"V_ (") = 0(n*) and m matrices of order 
Eto (?) = 0(n*-i). " 

2.3. A new block-diagonal hierarchy. One drawback of the Lasserre hierarchy 
is that the computational cost for optimizing over Q l (K) is considerably higher than 
the cost for optimizing over N^iK). To define a more economical variation of it, 
which still refines the Lovasz-Schrijver hierarchy, we consider a suitable principal 
submatrix of the full matrix Mt(y). 

For a positive integer t and a subset T C V of cardinality t — 1, let M(T;y) 
denote the principal submatrix of M t (y) whose rows and columns are indexed by 

A(T) := |J As, where As := {5*} U {S U {i} | i G V}. (4) 

SCT 
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It will be convenient to consider A(T) as a multiset: We keep possible repeated 
occurrences, e.g. S and S U {i} if i € S. So strictly speaking the matrix M(T;y) 
is a principal submatrix of M t (y) only after removing repeated rows and columns. 
We consider multisets here because it simplifies the notation in Lemma \2. 21 Note 
that M(0; y) = M\(y), and observe that we only need to know the components of 
y indexed by Vt+i{V), instead of T'2t(V) as in the Lasserre hierarchy, in order to 
define the matrices M(T; y) for all T G 7>=( t _i)(V). 

Define the first iterate of the block- diagonal hierarchy by L {K) := Q 1 ^) and, 
for t > 2, define its t-th iterate by 

L\K) := {x e K Pl(y) |3y e M p '+i( y ) : yo = z , w = Xi (i € V), 

M(T;;/)H(TeP =M (y)), 

M(T; o,y) h (T e V=(t-2)(V), 1 = 1,..., m)}. 

Since we used principal submatrices of the Lasserre hierarchy, we obviously have 
that the Lasserre hierarchy refines the block-diagonal hierarchy. As we see in Sec- 
tion 12.51 the block-diagonal hierarchy still refines the Lovasz-Schrijver hierarchy. 

Next we give a more compact formulation for the set based on the fact 

that the matrix M(T; y) has a special block structure which can be exploited to 
block-diagonalize it. This property justifies the name "block-diagonal hierarchy". 

For a subset S of T, let As(y) denote the principal submatrix of M(T; y) indexed 
by the set As, which is defined in <(5J). It is a (n + 1) x (n + 1) matrix lying in 
■M+ v with entries 

As(y)o,o = ys, A s (y) ,i = ysu{i}> A s (y)i,j = ysu{i,j} (*>ieV). 

The submatrix M(T; y) [S, S'} of M(T; y) with row indices in As and column indices 
in A s > depends only on5U S': M (T; y)[S, S'] = A Su s>(y)- 

Lemma 2.2. The matrix M(T;y) is positive semidefinite if and only if for all 
subsets S of T the matrix 

A(S,T)(y);= £ (-l)l s '\ s U s ,(t/) (5) 

S':SCS'CT 

is positive semidefinite. 

Proof. The proof is a "block-matrix version" of the one of {2} in [9|. Define the 
block-matrix Z indexed by A(T), whose (S, 5")-th block is the identity matrix I of 
order n + 1 if S C S' and the zero matrix otherwise. Its inverse is the block matrix 
whose (S, S")-th block is (— l)' s if S C S' and the zero matrix otherwise. Define 
the block diagonal matrix D with diagonal blocks A(S, T){y) for S C T. Direct 
verification shows that M(T; y) = ZDZ T . Therefore, 

M(T, y)hO ^> D>Q <S=^ V5CT: A(S, T)(y) h 0. □ 
Example 2.3. for T = {1, 2}, A(T) = A U A U A 2 U ^12 a«c? 



/A 


-4i 


A 2 


A 12 \ 




Ai 


A12 


A12 




A12 


A 2 


A12 


V^12 


A12 


A12 


A 12 J 



M(T;y) = 
where we wrote As instead of As{y). 



>- 



A - A x - A 2 + A 12 h 
At - A 12 t 
A 2 - A 12 y 
A12 h 
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Hence, in the formulation of L'(A), each condition M(T; y) >z0, which involves 
one matrix of order 2*~ 1 (n+l), can be replaced by the 2 t_1 conditions A(S, T)(y) >z 
0, each involving a matrix of order n + 1. Similarly, the condition M(T; aiy) >z 
can be replaced by the 2*~ 2 conditions A(S,T)(aiy) >z 0, each involving a matrix 
of order n + 1. 

2.4. A variation of the block-diagonal hierarchy. The next lemma deals with 
other possible ways of encoding the linear conditions defining the set A. It moti- 
vates our second variation L (A). It turns out that it has an explicit link to the 
Lovasz-Schrijver hierarchy. A main advantage of L l (K) over L'(A) is that we do 
not need an explicit linear description of the set A in order to be able to define 
L'(A). Hence £*(A) enjoys the same complexity property as N+(K): If one can 
optimize in polynomial time over A then the same holds for L (K) for any fixed t. 

Lemma 2.4. Let t > 1, y G M. Vt + l( - v \ A be as in flU and A(S,T)(y) be as m 0j. 
Then, the following two assertions are equivalent: 

(a) For all T G V ={t - 1} (V), S CT,ieV: 

A(S, T) (y) ei G A, A(S, T) (y) (e - e,) G A. 

(b) For all T G T =t (V), l = l,...,m: 

M T {aiy) h 0. 

Proof. Using the identities 

afA s (y)e = (aiy)s, aj A s {y)e l = (aiy) SU {i}, 
the conditions A(S,T)(y)ei G A, A(S,T)(y)(eo — e,) G A can be rewritten as 

J2 (-l) |SV W)s'u{i} >0 (l = l,...,m), 

S':SCS'CT 

]T (-V lS ' XSl ((aiy)s>~(aiy) S 'u{ t })>0 (i = l,...,m). 

S':SCS'CT 

On the other hand, using (J2]), Mr(aiy) >: is equivalent to 

£ (-l) |s '\ s W)s< >0 (SCT). 

S':SCS'CT 

From this one can verify the equivalence of (a) and (b). □ 

Observe that for t = 1 property (a) is equivalent to Ao(y)ei, ^4o(y)(eo — &i) G 
A for all i E V. Combined with the condition Ao(y) >z 0, this characterizes 
membership in the set N + (K). 

This motivates replacing in the definition of L*(A) the condition U M(T; aiy) >z 
for all T G 7>=( t -2)(V)" by property (a): For i > 1 define 

L\K) := {a; G R Pl(v) |3y G R^( v ) : y = xq,^ = a* (i G V), 
M(T;j/)bO (T6P =(t _i)(V)), 
A(S,T)(y) ei G A, A(5,T)(i/)(eo - e<) G A 
(T G V = ( t -i)(V), S C T, i e V")}. 
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2.5. Comparisons. Another advantage is that £'(AT) can be directly compared 
to the Lovasz-Schrijver hierarchy N\.{K). The next proposition shows that our 
second variation refines the Lovasz-Schrijver hierarchy. 

Proposition 2.5. We have L^K) = N+(K) and L\K) C N\{K) fort>2. 

Proof. As noted above we have L l (K) = N+(K). Now let t > 2 and x £ L l (K). 
Thus, there is a y £ K' p t+i( v ') which satisfies yo — xo, yi — Xi {i £ V), and 
M(T;y) h or, equivalently, A(S,T){y) h for all S C T C V with |T| = t - 1. 
Moreover property (a) of Lemma T2.4I holds. Set Y := Mi(y). Then a; = Yeo and 
Y £ .M-|-,y. Given 1 < s < t — 1, and ii, . . . ,i s £ V, and ct € {±1} S , consider 
the multisets T = {ij, . . . , 5 = {i r | r = 1, . . . , s, ay = 1} C T, and define 
Y^';;^ 8 := A(5,T)(y). Here we extend the definition of A(S,T)(y) in © to the 
case when S and T are multisets by taking the summation over all multisets S' 
lying between S and T; moreover, when 5" is a multiset with S" as underlying set, 
we let As'{y) := A$"(y). Now one can verify that the conditions from Proposition 
HO hold, which implies x € N\ (K). □ 

As one can see from the above proof, the main difference between i'(if) and 
N+(K) is that the matrices Y? 1 '"'!*' share many common entries in the definition 

of L*(i<r). As a consequence, one can describe the set L l (K) with less variables 
compared to N+(K). In Table 1 we compare the complexity of the formulations for 
I^iK) and N^{K). In both cases one has a semidefinite programming formulation 
involving a number of matrices of size n + 1 required to be positive semidefinite 
and a number of conditions of the type "x £ K" . 





L\K) 


iV|(Jf) 


# variables 


ESo (?) 


G)EU(2n) 1 




= 7^Ty T ^+ 1 + 0(n t ) 


= 2'- 2 n t+1 + 0(7i') 


# matrices 




E: =0 (2n) 1 


of order n + 1 


= & t " 1 + o(- t - 2 ) 


= 2 t " 1 n t - 1 + 0(n t - 2 ) 


$ conditions 






"x £ K" 







Table 1 . Complexity comparison for L l (K) and (K) . 



Also, as already stated in Section 12.31 the block-diagonal hierarchy refines the 
Lovasz-Schrijver hierarchy. This can be seen by comparing L t+1 (K) with the second 
variation L l (K). 

Proposition 2.6. For t> 1 the inclusion L t+1 (K) C L l {K) holds. 

Proof. This follows directly from the definitions, after noting that, for |T| = t, the 
index set of Mr(y) is contained in the index set of M(T\ {i};y), where i is any 
element of T. □ 

3. Application to the stable set problem 

In this section we apply the new hierarchies to the stable set problem. Let 
G = (V, E) be a graph. A subset S C V is called a stable set if none of its vertices 
are adjacent. The incidence vector of S is x £ {0; 1} with x S (*) = 1 iff i G S 1 . 
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The stability number a(G) is the maximum cardinality of a stable set. By Sg we 
denote the set of all stable sets of G. Then the stable set polytope is 

STAB(G) := conv{x S | S G S G }, 

and the corresponding cone is 

ST(G) :=R + {( X X S ) | SeS G }. 
A linear relaxation of ST(G) is the fractional stable set cone 

FR(G) := {x G K Pl(y) | a* > (i G V), x t + Xj < x ({i,j} G E)}. 
A semidefinite relaxation of ST(G) is the theta body 

TH(G) := {Ye G R Vl{V) \ Y G M+y,Y 3 = ({ij} G £)}, 

which is contained in FR(G). Maximizing the linear function X^ev x i over the 
theta body TH(G) intersected with the hyperplane xq = 1 equals the Lovasz theta 
function $(G) introduced by Lovasz in [12]. For details about these relaxations and 
the stable set problem we refer e.g. to [IT] and [15] . 

In [9] Lemma 20] it was shown that when constructing the Lasserre hierarchy for 
FR(G) one can considerably simplify the formulation. One can replace the condi- 
tion u M t -i(aiy) y 0", where ai runs through all linear inequalities defining FR(G), 
by the simpler equalities = 0, where {i,j} £ E, the so-called edge equalities. 
We want to apply the same simplification to the definition of L*(FR(G)) and de- 
fine another variant L (G) of it. However, in contrast to the Lasserre hierarchy, 
this simplification weakens the block-diagonal hierarchy a little bit since we can 
only claim the inclusion L*(FR(G)) C L*(G). Nevertheless the new variant £*(G) 
still refines the Lovasz-Schrijver hierarchy, as L l {G) C N t ^ 1 (TH(G)) follows from 
Proposition 13.21 below combined with Proposition ^. 51 We define 

L*(G) := {x G R Pl(y) \3y G R v ^ v *> : y = x , Vl = x. t (i G V), 
M(T;y)hO(TeV Ht . 1) (V)), 

yij = o({i,j}eE)}. 

Thus, L l (G) = TH(G) and one can easily verify the inclusions Z t (FR(G)) C L\G) 
when t > 1 and L*(FR(G)) C L*(G) when t > 2. Maximizing the objective function 
J2iev x i over L 2 (G) intersected with the hyperplane xo = 1 coincides with the 
parameter 1(G) considered in lol [71 ITU], 

The next lemma says that the edge conditions in the definition of £*(G) imply 
that all variables indexed by non-stable sets are identically 0. 

Lemma 3.1. Let y G M.' Pt+1 ^ v ' satisfy the conditions in the definition o/L*(G). 
Then yi = for any subset I C V with \I\ < t + 1 and containing an edge. 

Proof. For |/| = 2 the statement is nothing else but the edge equalities. Assume 
that | J | > 3, let i, j G I be adjacent vertices, and let k be another vertex in /. Define 
T := I \ {i,j}- The matrix M(T; y) is positive semidefinite and the sets {i,j} and 
T U {k} occur in the index set A(T). As the (ij, ij)-th entry of M(T; y) is — 0, 
we have by the positive semidefiniteness of M(T;y) that its (ij, T Li {fc})-th entry 
is as well and the statement of the lemma follows. □ 

Proposition 3.2. We have the inclusion L t+1 (G) C Z*(TH(G)) fort>\. 
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Proof. Assume that y G R' Pt + 2 ( v ) satisfies the conditions of the definition of L t+1 (G). 
In the following we show that the vector consisting of the first n + l coordinates of 
y belongs to Z*(TH(G)). 

Fix T € V t -i(V), S C T and k e V. We show that A(S,T)(y)e k £ TH(G) and 
A(S,T)(y)(e - e k ) € TH(G). For this we construct matrices Y k and Z k in 
satisfying Y% = Z k j = when i and j are adjacent, and satisfying 

Y k e = A(S,T)(y)e k , Z k e = A(S,T)(y)(e - e k ). 
We distinguish between three cases. 

(1) k £ S: Then A(S,T)(y)e = A{S,T)(y)e k and define Y k := A(S,T)(y), 
Z k := 0. 

(2) keT\S: Then A{S,T)(y)e k = and define Y k := 0, Z k := A(S,T)(y). 

(3) keV\T: Then we define Y k = A(S U {k},T U {k})(y) and Z k : A(S,TU 
{k})(y) = A(S,T)(y) - A(SU{k},TU {k})(y). 

In all cases we see by Lemma [2~2l 13 . 1 1 that Y k , Z k satisfy the desired conditions. □ 
We summarize the inclusion relations between the various relaxations: 
ST(G) C Q'(FR(G)) C i*(G) C L^iTUiG)) C Af^T 1 (TH(G)). 
Moreover, iV+ _1 (TH(G)) = ST(G) holds for i > a(G) (see for a proof). 

4. Experimental results 

In this section we present some computational results for Paley graphs. 

Let ¥ q be the finite field with prime power q which is congruent to 1 modulo 
4; then —1 is a square in ¥ q . The Paley graph P q has ¥ q as vertex set and two 
distinct elements u, v € F 9 are adjacent if u — v is a square in ¥ q . The Paley 
graph is isomorphic to its complementary graph, it is a strongly regular graph and 
its automorphism group acts doubly-transitive on the vertices. It is known ( |12[ 
Theorem 8]) that i?(G)i9(G) = |V(G)| when G is a vertex-transitive graph, where 
G denotes the complementary graph of G. Since the Paley graph P q is vertex- 
transitive and isomorphic to its complementary graph, we have $(P q ) = $(Pq) = 
y/q (cf. [U Theorem 13.14]). J.B. Shearer ([TB]) has computed ct(P q ) for all primes 
q < 7000. For more information about P q we refer e.g. to [2j Chapter 13.2]. 

In order to illustrate the quality of the new relaxations L 1 (Pq), we have computed 
the bounds obtained by maximizing X)ugy(p ) x v over the sets Z/'(P g ) (for t = 2, 3) 
and N + (TH(P q )) intersected with xq = 1. The results are given in Table 2. There 
we consider all primes q congruent to 1 modulo 4 between 61 and 337, as well as a 
few larger values of q up to 809. We have chosen the Payley graph here because its 
automorphism group acts doubly-transitive on the vertex set and so our formulation 
for L\P q ) (t < 3) and N + (TR(P q )) considerably simplifies. (See Chapter 6.1] for 
implementation details.) For instance, optimization over L 3 (Psog) (resp., L 2 (P%q$), 
A r + (TH(P 80 9))) can be formulated via an SDP with 876 (resp., 36, 812) variables 
and with four matrices with sizes 808, 808, 404 and 202 (resp., two matrices with 
sizes 809 and 405, three matrices with sizes 810, 810 and 809). For the computations 
we used the program CSDP [3|. Experiments were conducted on a single machine 
with an Intel(R) Pentium(R) processor, 3Ghz and 1GB of RAM. To compute the 
bounds from Table 2 we needed less than a minute when q < 100 and, for the 
largest instance Pso9, around 45 minutes for L 3 (P 80 g), 31 minutes for L 2 (P 8 og) and 
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4.5 hours for iV+(TH(P809))- Thus as expected the relaxation L 2 (G) gives a sharper 
bound than _/V+(TH(G)), however at a much smaller computational cost. 

Finally note that one can strengthen the relaxation L f (G) by adding the non- 
negativity constraints y > 0. However this only gives a marginal improvement for 
Paley graphs, as the bounds differ only in decimals. 



1 


L 1 (P q ) = TU(P q ) 
■d(P q ) = 


Ar+(TH(P 9 )) 


L 2 (P q ) 




a(P q ) 


61 


7.810 


5.901 


5.465 


5.035 


5 


73 


8.544 


6.377 


5.973 


5.132 


5 


89 


9.434 


7.155 


6.304 


5.391 


5 


97 


9.849 


7.948 


7.398 


6.596 


6 


101 


10.050 


7.290 


6.611 


5.496 


5 


109 


10.440 


8.007 


7.366 


6.578 


6 


113 


10.630 


8.330 


7.599 


7.009 


7 


137 


11.705 


8.829 


8.200 


7.047 


7 


149 


12.207 


9.188 


8.231 


7.136 


7 


157 


12.530 


9.695 


8.707 


7.485 


7 


173 


13.153 


10.316 


9.426 


8.062 


8 


181 


13.454 


10.324 


9.112 


7.606 


7 


193 


13.892 


10.506 


9.210 


7.651 


7 


197 


14.036 


10.652 


9.226 


8.064 


8 


229 


15.133 


11.659 


10.290 


9.076 


9 


233 


15.264 


12.382 


10.182 


8.245 


7 


241 


15.524 


11.595 


9.891 


8.275 


7 


257 


16.031 


11.558 


10.247 


8.131 


7 


269 


16.401 


12.307 


10.624 


8.778 


8 


277 


16.643 


12.469 


10.340 


8.670 


8 


281 


16.763 


11.902 


10.605 


8.397 


7 


293 


17.117 


13.127 


10.937 


9.183 


8 


313 


17.692 


13.128 


11.630 


9.458 


8 


317 


17.804 


13.861 


12.377 


10.375 


9 


337 


18.358 


13.724 


11.658 


9.464 


9 


401 


20.025 


14.927 


12.753 


10.023 


9 


509 


22.561 


16.580 


14.307 


11.196 


9 


601 


24.515 


17.999 


16.077 


12.484 


11 


701 


26.476 


19.332 


16.857 


12.822 


10 


809 


28.443 


20.636 


17.371 


13.499 
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Table 2. Optimizing over L*(P q ) and N+(TK(P q )) for Paley graphs. 
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